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Abstract 

We study the critical properties of the spin-continuous S 4 system on the typical translational in- 
variant triangular lattices by combining the recently-developed generalized Migdal-Kadanoff bond- 
moving recursion procedures with the cumulative expansion technique. In three different cases of 
nearest-neighbor, next nearest neighbor and external field we obtain the critical points and further 
calculate the critical exponents according to the scaling theory. In all case it is found that there 
exists three fixed points. The correlation length critical exponents obtained near the Wilson-Fisher 
fixed points are found getting smaller and smaller with the increasing of the system complexity. 
Others are found similar to the results of the classical Gaussian model and different from those of 
the Ising system. 
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I. INTRODUCTION 



The S 4 model, which is another extension of the classical Ising model, is a kind of typical 
spin-continuous systems which allows the spin to take any real value between (—00, +00) 
instead of a discrete one. The study of it is of great value for better understanding the 
properties of the ferromagnetic systems from either a theoretical or a practical viewpoint. 
Due to this point its it has been widely studied in the past few years [l-3]. For an example, 
critical properties of the S 4 model on some lattice systems with translational symmetry has 
been found to have close dependent on the space dimensionality d : when d is more than 4, 
only the Gaussian fixed point is obtained while for low-dimensional systems with d less equal 
to 4, not only the Wilson-Fisher fixed point but also the Gaussian one can be obtained 
Similar dimensionality dependent behaviors are also found for the S 4 model on the fractal 
lattices 

However, critical behaviors of the S 4 system is far from having been well studied. This 
may in some case due to the self- complexity of the model caused by its four-body interac- 
tions. Other restrictions may come from the particular lattice system on which the model 
is constructed. Recently, a type of generalization of the remarkable Migdal-Kadanoff bond- 
moving renormalization procedure is developed which enables us to handle complex systems 
in a relatively easy way (4J. For the great convenience it has brought in the coarse-graining 
process, in this paper, we extend this renormalization procedure to study the critical prop- 
erties of the S 4 system on translational invariant triangular lattice by combining it with the 
cumulative expansion method {2, 3, 5|. 

The paper is organized as follows: in Sec. [TTJ a brief review of the generalized bond- 
moving recursion procedures is presented; in Sec. IHIt the procedures are used to study the 
critical properties of the spin-continuous S 4 model constructed on the translational invariant 
triangular lattice; Sec. IIVI serves as a summary of our conclusion in which further implicit 
applications are also discussed. 

II. RENORMALIZATION PROCEDURE AND CUMULATIVE EXPANSION 
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In Ref. |4|, a new type of generalization of the remarkable Migdal-Kadanoff bond- moving 
renormalization procedure is developed bringing with much convenience for the study of 
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FIG. 1: Generalized bond- moving procedure recurring on the triangular lattice where the peripheral 
bonds connecting the three to be eliminated sites 1, 2 and 3 around the selected triplet A ABC are 
drawn at half length. 

spin-continues systems. By introducing in some kind of symmetrical half-length bond oper- 
ations, the renormalization process is greatly simplified as has been witnessed in the study 
of Gaussian system. For the example of systems constructed on the triangular lattice as il- 
lustrated in FigJU this was realized by supposing each bond connecting the to be eliminated 
sites denoted as 1, 2 and 3 in A ABC to schlep symmetrically two half-length bonds. Thus 
the number of lattice sites that a basic unit for recursion contains can be reduced to only 
six which provides a good foundation for the following integral calculation. 

The technique of cumulative expansion is always used as an auxiliary method in calculat- 
ing high order integrations {2, 3, 5|. For a system with effective Hamiltonian: H c g = Hq + V, 
a partial trace (PT) is defined as 

( pT ) = r (n^) eHefl 
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is namely the cumulative expansion average and 

p+oo f m 



The integrations in the above equations trace over all the to be decimated sites in the 
renormalization process. 

On the supposition of V is a small variable, we can expand e v by a Maclaurin series as 

e v = l + V+±V 2 + ±V 3 + ---. (4) 

The partial trace is then obtained to be 

(PT) =A(l + (V)o + ^(V\ + i(y 3 ) + ••-). (5) 

Noticing that the partition function should be kept unchanged after the renormalization 
group (RG) transformation. That is, supposing 



/+oo / m \ 
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(6) 



then H' eS represents an effective Hamiltonian after RG transformation (C is a RG constant). 
Thus we can obtain from the above equations 

H' eS = \nA + (V) + \ {(V\ - (V)l) + ■ • • , (7) 

2 3 

where the approximation relation ln(l + x) = x — \ + \ — • ■ • is used. 

Now we have found a convenient way to calculate the effective Hamiltonian H' eS from the 
partial trace defined herein before. Supposing 
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JI^ )e H ° = e H °, (8) 



the first part in the right-hand site of Eq.Q can be regarded as the zero order approxi- 
mation of H' eS . Similar understandings are still true as well for other right-hand parts. In 
the following sections, this will be applied to study the critical behavior of the S 4 system 
constructed on the triangular lattice combining with the generalized bond-moving renormal- 
ization procedures. 
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III. CRITICAL BEHAVIOR OF THE S 4 SYSTEM 



The S model is another spin-continuous extension of the Ising model [6| . By introducing 
a probability distribution function 

W( Sl ,s 2 ,--- ,s N ) = e- b ^ s i-^ us ", (9) 

for the classical Ising spin system, the effective Hamiltonian of the S 4 model in an external 
field h can be obtained as 

H eS = Y,Ks iSj -~Y,Si-uY,4 + hY, s " ( 10 ) 

(ij) i i i 

where K = J / ksT is the reduced interaction with K > denotes the ferromagnetic system; 
ks is the Boltzmann constant and T the thermodynamic temperature; h is the external field. 
The (ij) in the first summation generally represents a certain nearest-neighbor spin pair. 
Seen from Eq. ffTO]) . the S 4 model Hamiltonian contains synchronously the four-spins and 
two-spins interactions denoted by u (u > 0) and Gaussian constant b respectively. This will 
undoubtedly made it more difficult to be studied, however these can all be settled following 
the generalized bond-moving renormalization procedures. 

A. nearest-neighbor interaction 

Firstly we give our study on the case of only considering nearest-neighbor interactions. 
The procedure of the generalized bond-moving RG transformation for obtaining a certain 
bond K' is shown in Figj2] as an example. From which we can see that the number of to be 
decimated sites reduces to be only one. Thus the effective Hamiltonian of the basic unit for 
RG transformation in the case of h = can be given as 

H cS = 2K {s aSl + Sl s b ) - b (2s\ + s 2 a + si) - 2u (2s 4 + s 4 + s 4 b ) , (11) 

where in the RG procedures two types of interactions K e and K together with two types of 
self-energy (—b e s 2 /2) and (—bs 2 /2) are assigned respectively for differentiation of the to be 
and not to be eliminated bonds Gefen et al did on the fractal 7H9] . For the particular case 
of triangular lattices the numerical value of K e and K is actually identical as well as that of 
b e and b (or u e and u). In the recursion procedures the two half length bonds are considered 
acting effectively as a whole one and be moved regularly as other bonds. 
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FIG. 2: Bond-moving and decimation procedures for the renormalized bond K' between sites A 
and B on the triangular lattices. 



Supposing 



and 



V = -lus\, 



H = 2K (s aSl + s lSb ) - b (2s\ + s 2 a + si) - 2u (s* + sf) 



(12) 
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the zero-order approximation of the effective Hamiltonian after RG transformation can be 
obtained from Eq.flS]) as 



H' = K 0l s a s b + K 02 (s 2 a + si) + K m (s 4 a + 4) , 



with 
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Continuously we can obtain the first-order one by Eq.([2|) as 



= K llSa s b + K 12 (s 2 a + s 2 b ) + K l3 (st + si) 



(16) 



where 
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In the same way, the second-order term can also be approximately obtained as 
\ ((^ 2 )o - (V)l) = K 21 s a s b + K 22 (s 2 a + si) + K 23 (s A a + st) , 



with 
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Combining Eqs. ffTo"]) . f[T7|) and 
mat ion 



we obtain the effective Hamiltonian after RG transfor- 



H' eS = (K Q1 + K n + K 21 )s a s b + (K 02 + K l2 + K 22 ) (s 2 a + s 2 b ) 
+ (K 03 + K 13 + K 23 ) (si + si) , 



(20) 



In order to get the recursion relations of the RG transformation, we rescale the spins by 
s i = £i s i = a -> b) and then Eq. (f2~TI) can be rewritten as 



K S = K's'J b - - (s' 2 + s' 2 ) - u> « 4 + sf) 



where 
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Setting K' = K = K* and u' = u = u*, the fixed points of the system can be found are 



A : K* = 0, u* = 0; (23a) 
B : K* = 6, w* = 0; (23b) 
C : K* = 0.4366, u* = 0.3706 2 ; (23c) 

in which A is found to be a nonphysical steady fixed point while B is a Gaussian fixed point 
and C a Wilson-Fisher one. 

Actually, in the three fixed point, the Wilson-Fisher one is most crucial for the critical 
properties of the S 4 system. By expanding K' and u' in the neighborhood of the Wilson- 
Fisher fixed point C and preserving only the linear items, the transformation matrix is 
obtained to be 

/ % mi \ ( 2.855 2.895 \ 
Rl(K,u)= = , (24) 

\lk d ik) c v 1 - 356 2M7 J 

with two eigenvalues Ai = 4.643 > 1 and A2 = 0.659 < 1 revealing that C is a barrier 
point. Thus a scale power can be obtained according to the scaling theory p = = 1.107 
where d = 2 is the dimensionality of the triangular lattice and L = 2 is the scaling factor. 
The correlation length exponent can be obtained from the relations u = ^ to be y = 0.451 



smaller than the Gaussian case which contains only two-body interactions [4j, [10] while in 
good conformity with previous S 4 system studies on other lattices as. 



B. external field 

In order to obtain all the critical exponents, we study in the following the case of S* 4 
system in an external field h where the effective Hamiltonian can be written as Eq.f JTQ|) . For 
simplicity, we give out here only the primary results while leaving the process of derivation 
in the appendixes. Considering only nearest-neighbor interactions, three fixed points can 
also be obtained by following similar RG transformation procedures. They are 

A: K* = 0,u* = 0, h* = 0; (25a) 
B: K* = b, u* = 0, h* = 0; (25b) 
C : K* = 0.4366, u* = 0.35 26 2 , h* = 0; (25c) 
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the transformation matrix in the neighborhood of the Wilson-Fisher fixed point C is now 

/ o C70 o n90 n \ 



( dW_ dIC diC_ \ 

OK du dh 

du' du' du' 

dK du dh 

dh' dh' dh' 

\ dK du dh / 



C 



2.872 3.032 
1.307 2.450 
2.875 



(26) 



with three eigenvalues Ai = 4.663 > 1, A 2 = 0.659 < 1 and A 3 = 2.875 > 1. Where Ai and 
A2 is related to the temperature while A3 is related to the external field. Then we can obtain 
the two scale powers according to the scaling theory 

InAi hiA 3 



P 



d\nL 



1.110, 



(ML 



0.762: 



The critical exponents at this point can be obtained following a = , (5 = 7 



(27) 

2g-l 



5 = j^, r] = 2 + d(l - 2q), v = f d to be a = 1.099, (3 = 0.215, 7 = 0.471, 5 = 3.197, 
7] = 0.953, v = 0.450 describing in the neighborhood of the critical point the variations 
in heat, magnetization, magnetic field, magnetic susceptibility, correlation function and 
correlation length, respectively. 



C. next-nearest neighbor interaction 

For the enriching of the S A results, we present in this paper also the case of considering 
next-nearest neighbor interactions. All the approximation of the effective Hamiltonian after 
RG transformation can be obtained also by combining the generalized bond-moving opera- 
tions and the cumulative expansion techniques. The only change need to be made is to select 
a cluster of containing several sites (as is shown in Figj3]) as a basic unit for recursion. The 
various coefficients of H' eS after RG derivation is also left in the appendixes for simplicity, 
we give out here only the primary results. 

In the case of considering next-nearest neighbor interactions, we can also obtain by RG 
transformation three fixed points 

A: K* = 0, u* = 0; (28a) 
B : K* = 0.8296, u* = 0; (28b) 
C : K* = 0.2576, u* = 5.4056^; (28c) 

in which A is found to be a nonphysical steady fixed point while B is a Gaussian fixed point 
and C a Wilson-Fisher one. After some algebra we can obtain the correlation length critical 
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FIG. 3: Bond-moving and decimation procedures for the renormalized bond K' in the case of 
considering next nearest-neighbor interactions. 

exponent at the Wilson-Fisher fixed point to be v = 0.341 also in conformity with previous 
studies. However we find it becomes relatively smaller than previous results. This reveals 
the effects of next nearest neighbor interactions. 

IV. SUMMARY AND DISCUSSIONS 

In summary, we have presented in this paper a detailed study on the critical behavior 
of the S 4 system constructed on the triangular lattices by combining the generalized bond- 
moving RG transformation method and the cumulative expansion techniques. In all cases it 
is found to have three fixed points. The critical exponents obtained near the Wilson-Fisher 
fixed points are found relatively smaller than that of the Gaussian model and different from 
those of the Ising model. This reveals the particular effect of the four-body interactions on 
the critical properties of the S* 4 system. In further by comparing all the results of correlation 
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length exponent obtained in all cases we find that it becomes smaller and smaller with the 
increasing of the complexity of the system. This reveals the decisive influence of the system 
complexity on the critical behavior. 
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K M = (A4) 



Appendix A: RG coefficients of H' eff in external field 

For simplicity we give out here only the various RG coefficients of H' cS from zero to second 
order approximations instead of showing all the details of derivation that is very similar to 
those have appeared in this paper hereinbefore. In the case of considering the influence of 
an external field h, these coefficients are 

K 2 

Koi = — , (Al) 

*. - ^ m 

K 03 = -2u, (A3) 
K -2b 
b~ 

Kn = - 3J ^(b + h 2 ), (A5) 

Kn = -^{b + h% (A6) 

Kn = (A7) 
Kvh 

Ku = -^{Zb + h% (A8) 
3iCV /or 93/i 2 9h 4 h\ 

K » = ^ 35 + — + 1F + ^ (A9) 

3K 2 u\ 93h 2 9h 4 h\ , k . 

^ = - Si -(86 + — + + (aw) 

A' 4 u 2 „ 174ft 2 17ft 4 , ,,,,, 

K » = w (87+ — + — >■ (A11) 

hKu 2 „ nr 105/i 2 21h 4 h\ . Ai . 

K ™ = "1^ (105 + — + ~W + ¥ y (A12) 
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where K i4: (i = 0, 1,2) represents the various RG coefficients of H' cS which is related with 
the external field. 



Appendix B: RG coefficients of next-nearest H' eff 

The RG coefficients for various order approximation of the effective Hamiltonian in the 
case of considering next-nearest neighbor interactions are 

48K 2 b + 51K 3 

Km ~ 16(462 _ KiT) ' ™ 

966 3 - 75K 2 b - 12K 2 
K ° 2 ~ 16(4^ — k*) ' ^ 

^os = -3m, (B3) 

Ku ~ 463(462-^)2' 

K 2 u{272b 4 + 128kb 3 + 136K 2 b 2 - 17K 4 ) 
Kl2 ~ 6463(462 - K 2 ) 2 ' (B5) 

K 4 u(±\\2b 4 + 217QKb 3 + 2824K 2 6 2 + 1088K 3 6 + 771K 4 ) 
13 ~ 40 9 66 4 (46 2 - K 2 ) 2 ' ( ' 

2K 2 u 2 (Q5Qb 5 + 180 26 4 fT + 452b 3 K 2 + 374^3 + mK ± + llK ^ 
K * ~ 9(462 -^)5 > v G 

fr 2 w 2 (25846 5 + 14086 4 ir + 12586 3 fT 2 + mb 2 K 3 + ll%K 4 + 8K 5 ) 
K22 ~ 9(462 - K 2 f ' (B8) 

K 4 v 2 

K23 = 24 5 76(463 - bK 2 f (157900806 12 + 94699526^ + 60375046^ 2 
+54312966 9 fsT 3 + 44208646 8 J fsT 4 + UUU2b 7 K 5 + 1850886 6 fT 6 
-20 8 8 9 66 5 ^ 7 - 1525926 4 fT 8 + 130 5 66 3 fsT 9 + 208086 2 ^ 10 - 8Q7K 12 ). (B9) 

where the next nearest-neighbor interaction is set one quarter in strength as the nearest- 
neighbor interaction. 
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